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Abstract: We study asymptotically slowly varying perturbations of the AdS black brane in 
! Einstein's gravity with a negative cosmological constant. We allow both the induced metric 
Q ! and the Brown- York stress tensor at a given radial cut-off slice to fluctuate. These fluctuations, 
^ ! which determine the radial evolution of the metric, are parametrized in terms of boundary data. 
'. We observe that the renormalized energy-momentum tensor at any radial slice takes the stan- 
! dard hydrodynamic form which is relativistically covariant with respect to the induced metric. 
2 ! The RG flow of the fluid takes the form of field redefinitions of the boundary hydrodynamic 
^ i variables. To show this, up to first order in the derivative expansion, we only need to investigate 
^ I the radial flow of the boundary data and do not need to impose constraints on them. Imposing 
^ I the constraints gives unforced nonlinear hydrodynamic equations at any radial slice. Along 
I the way we make a careful study of the choice of counter-terms and hypersurfaces involved in 
defining the holographic RG flow, while at the same time we do not explicitly set any boundary 
condition either at the cut-off or at the horizon. We find that t]/s does not change along the 
RG flow, equaling l/(47r) when the future horizon is regular. We also analyze the flow of the 
speed of sound and find that it diverges at the horizon. 
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1. Introduction 

Since the early days of AdS/CFT it has been intuitively clear, that the radial direction 

which takes us away from the boundary of the space-times dual to the CFT states, encodes 
the CFT's renormalization group (RG) flow. However, a clear and general formulation of how 
the evolution of the boundary variables in the radial direction as given by Einstein's equations, 
realizes the field-theoretic renormalization group fiow, still eludes us. 

Instead of trying to gain a full understanding of the holographic renormalization group fiow, 
here we will try to test the correspondence in a regime where the gravitational fiow equations 
become more tractable. 

Such a simplification indeed occurs in the long-wavelength regime, namely by considering 
only low frequency fiuctuations in the field theory directions. The Einstein equations then can 
be solved order by order in the derivative expansion, and the results may be compared with the 
hydrodynamical expansion in the dual field theory. This framework is known as the fluid/gravity 
correspondence. 

One might argue, though, that it is a bit counter-intuitive to talk of RG fiow in the hydro- 
dynamic regime, since it is already a low-energy approximation. However, there is a natural 
description of a renormalization fiow of a non-linear partial differential equation like the Navier- 
Stokes' equation. It is known as Barenblatt's renormalization group [|, § and it's basic definition 
is as follows. Suppose a field variable follows a non-linear partial differential equation. When 
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u^'(x) = I d^x'u^ix') . (1.1) 



the field is coarse-grained (averaged) at a given spatial and temporal resolution, this averaged 
field follows a similar partial differential equation but with different parameters described by an 
appropriate renormalization group flow. Constructing this renormalization group flow is math- 
ematically a diflicult problem and can be easily implemented only for Hamiltonian systems 
nevertheless, numerical evidence suggests that it can be implemented for KdV-Burgers' 
and non-relativistic Navier-Stokes' equations as well ^ . This renormalization group flow gives 
key insights into turbulence. 

In the context of relativistic Navier-Stokes' equation, we are not aware of any attempt for 
renormalization maintaining manifest relativistic covariance. However, let us argue that in order 
to implement this we also need to evolve the space-time metric along the RG flow. Suppose, 
we average the four-velocity vector fleld over a space-time cell Vx around a space-time point 
X as below: 

1 

The averaged fleld u^{x) will not satisfy the normalization condition u^h^yu" = — 1 with h^j^y 
being the original metric in which the fluid is living. In order to satisfy this we need to 
consistently deflne an averaged metric as well. The holographic renormalization group flow 
achieves this beautifully as it naturally integrates the scale dependent hydrodynamic variables 
and the background d-dimensional metric together into the metric of a -|- 1 space-time. It 
will be interesting to explore in the future if the holographic RG indeed gives a realization of 
Barenblatt renormalization group flow of the relativistic Navier-Stokes' equation. However, here 
we will not attempt to make a fleld theoretic interpretation of the holographic renormalization 
of the fluid beyond these observations. 

In this paper we study a special class of geometries which are asymptotically both AdS 
and slowly varying, and tend to a boosted black brane solution far in the past or in the 
future. In the literature, such space-times have been studied by promoting the boost four- 
vector, the Hawking temperature and the flat metric on radial slices characterizing these black 
brane solutions, to arbitrary slowly varying functions of the boundary coordinates. Using the 
Dirichlet boundary condition at the boundary and regularity at the future horizon, one gets a 
systematic fluid/gravity correspondence H, ||, |T0|, [Tl|, |T2|, 0, |1^, Here we allow both the 



induced metric and the Brown- York stress tensor at a chosen radial slice to fluctuate and flnd 
how we can capture these fluctuations in terms of the boundary data. This implies that we are 
not imposing any explicit boundary condition on the radial flow of the boundary data. 

In accordance with expectations, the radius of a chosen radial slice is interpreted as the 
scale at which the observer probes the dual fleld theory. The product of this chosen radius po 
and the temperature fleld T{x), i.e. poT{x), is a conformal invariant. The data at p = po, which 
completely determines the radial flow, can be encoded at the boundary in terms functions of 
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this invariant. The shding energy scale of the renormahzation group flow becomes the variable 
cut-off radius that throughout this paper will be denoted by p^. 

We study the renormalized Brown- York stress tensor of these space-times at an arbitrary 
cut-off p = p'^ and flnd that, up to flrst order in the derivative expansion, we can always flnd 
field redefinitions of the hydrodynamic variables which brings the energy-momentum tensor 
to the standard (Landau-Lifshitz) hydrodynamic form, without using the constraints imposed 
by Einstein's equations on the data at any radial slice. We can find these field redefinitions 
by solving the radial fiow of the boundary data only. Thus the renormahzation group fiow of 
the fiuid takes the form of field redefinitions of hydrodynamic variables. Also, imposing the 
constraints of Einstein's equations implies the usual unforced nonlinear hydrodynamic equations 
at any cut-off radial slice p^. Further rj/s is the same at any p*^ irrespective of the explicit 
choice of boundary conditions and equals =Fl/(47r) when the past/future horizon is regular^. 
We also find that the speed of sound diverges at the horizon just as in the membrane paradigm 



0, o, 0, ig, 0, 

For the sake of explicitness, we will consider a five-dimensional space-time dual to a four- 
dimensional fiuid throughout this paper, though our results can be readily generalized to any 
dimension. 



In recent literature [^, ^ ^ (see also |32| for study of scale dependence of 
holographic response functions), the so-called "Wilsonian approach to fiuid/gravity duality" 
has already been investigated. However, some or all of the following features recur, which 
are assumptions of (z) the Dirichlet boundary condition at the cut-off as in the traditional 
holographic renormahzation, (ii) regularity at the future horizon as in AdS/CFT at finite 
temperature, (Hi) specific algebraic properties of the space-time,^ (iv) a non-relativistic limit, 
and {v) the linearized approximation. We show here these above-mentioned assumptions are not 
necessary in getting the essentials of the fiuid/gravity correspondence at any scale though they 
could be necessitated by other considerations. Moreover, rj/s is determined by the regularity 
at the past /future horizon alone at any scale irrespective of the explicit boundary condition at 
the cut-off. Our study is however restricted to asymptotically AdS space-times only. 

The motivation to remove the restriction of the Dirichlet boundary condition comes from 
some other recent investigations on the holographic renormahzation group flow in the literature. 



^We also point out that it is not completely irrelevant to talk about boundary conditions where fluctuations 
about equilibrium die in the past rather than in the future. In real time finite temperature field theory (i.e. 
in the Schwinger-Keldysh formalism) for instance, for local diagrammatic rules to work we need to take into 
account both the retarded and advanced propagator (along with the Keldysh propagator). The hydrodynamic 
limit of the advanced propagator has a pole in the upper-half plane corresponding to the flip in the sign of 
the viscosity. In real time gauge/gravity duality (which is still not a fully developed subject) one may expect 
similarly that the time-reversed boundary condition where the past horizons are regular should play a role. 

^In , a Petrov type 1 condition has been imposed on the geometry to reduce the bulk degrees of freedom 
of gravity to the fluid variables at the cut-off. 
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In particular, it has been argued that the counter-term action at the cut-off hypersurface in the 
traditional holographic renormalization prescription [0, |3^, |3^, ^ should be replaced 



by a boundary action which is the result of integrating out the bulk degrees of freedom in 
the region between the boundary and the cut-off ^ (see also In the saddle point 

approximation this boundary action should be obtained in conjunction with the bulk solution 
ro| and should not be just derived from the latter. The boundary action thus may affect 



the boundary condition at the cut-off. Making comparison with field theoretic renormalization 
scheme, a precise prescription can be given for a bulk free scalar field theory. 

However, it is not clear to us how to adapt these approaches to gravity. Nevertheless, we 
will do a careful study of the choice of counter-terms and hypersurfaces involved in defining the 
holographic renormalization group flow. We will show that in the context of the hydrodynamic 
energy-momentum tensor, by making some simple and reasonable assumptions, we can fix the 
choice of the counter-term and the hypersurface for the renormalization group flow of the holo- 
graphic fluid at the flrst order in the derivative expansion. In our approach we parametrize 
arbitrary boundary conditions at the cut-off, but it could be that the counter-term implicitly 
necessitates a deflnite boundary condition at the cut-off. We will investigate this issue fur- 
ther in future research, however for present purposes, the counter-term suffices to define the 
renormalized energy-momentum tensor at the cut-off. 

In our approach, we also do not need the non-relativistic limit of the bulk solution. Our 
results are however, consistent with the non-relativistic limit, as such a limit of any relativisti- 



cally covariant fluid gives incompressible Navier-Stokes equation [42|. This has been obtained 



holographically in the literature at an arbitrary radial slice as well [^, pq| . 

We also drop the requirement of regularity at the horizon, in order to work with arbitrary 
boundary data. In the context of AdS/CFT, we should restore the regularity condition at the 
end of our calculation, so that the transport coefficients get uniquely specified. Our general 



approach may shed light on the possible holographic emergence of gravity itself 44, 45, 46 



by relating holographic renormalization group fiow to absence of naked singularities in the bulk 
space-time. 

Finally, it is to be noted that the choice of bulk coordinates plays an important role in 
defining the renormalization group fiow. This is because the Brown- York stress tensor on a 
radial slice which gives us the scale dependent thermodynamic quantities and hydrodynamic 
transport coefficients is not covariant under transformation of the five bulk coordinates. So, 
choice of a bulk coordinate system is indeed a choice of a holographic renormalization scheme. 
In this paper, we choose the Fefferman- Graham (FG) coordinates. This coordinate system 
can be universally adapted for any asymptotically AdS space-time. Also, the one-to-one cor- 
respondence between conformal transformations in the boundary and bulk diffeomorphisms is 
obtained in any such space-time by finding the bulk diffeomorphisms which preserve the FG 
form of the metric. Therefore the radial coordinate in the FG coordinate system is a good 
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candidate for the field theory energy scale. In the specific cases relevant for the fluid/gravity 
correspondence, the bulk metric in the FG coordinates has a coordinate singularity at the 
horizon where the effective temperature diverges. This coordinate singularity signals an end 
of the renormalization group flow and may have an important physical significance. This has 
already been pointed out in where an attempt has been made to draw a parallel between 
holographic renormalization and entanglement renormalization by pointing out that both ends 
at a finite scale typically where the effective temperature diverges. In the latter case, it signifies 
that at this scale various parts of the system are disentangled from each other and we do not 
need to coarse grain any further. 

It should be stressed, that at any finite distance form the horizon the derivative expansion 
of the metric has a non-zero radius of convergence, and so we can trust the analysis in the FG 
coordinates, despite the coordinate singularity on the horizon. We also check the regularity 
of our gravity solution by transforming it to the Eddington-Finkelstein coordinates, which are 
regular at the horizon. 

Organization of the paper 

In section 2, we outline the asymptotically slowly varying AdS space-times and how they encode 
perturbations of both the induced metric and the Brown- York stress tensor at the cut-off radial 
slice. In section 3, we investigate the choice of counter-terms and hypersurfaces involved in 
defining the holographic renormalization group flow. In section 4 we find the field redefinitions 
which bring the renormalized energy-momentum tensor to the standard hydrodynamic form 
and also find the transport coefficients at an arbitrary radial slice. Finally, we discuss some 
open questions which we will like to investigate in the future. In the appendix, we present the 
gravity solution in the Eddington-Finkelstein coordinates. 

2. The asymptotically slowly varying AdS space-times 

We will be studying space-times which are: 

• asymptotically both AdS^ and slowly varying, 

• tend to the boosted black brane solution far in the past or in the future. 

Such space-times in Einstein's theory of gravity can be characterized by boundary data, which 
are the boundary metric and the boundary energy-momentum tensor. More specifically, the 
boundary energy-momentum tensor is the Balasubramanian-Krauss stress tensor, which is the 
appropriately renormalized Brown- York stress tensor when the cut-off is taken to the boundary. 
We will denote the boundary metric as g^^, and the boundary energy-momentum tensor as t^^,. 
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Instead of characterizing the space-times by boundary data, we will choose to characterize 
them by the data at a radial slice po- These data will be the induced metric and the Brown- York 
stress tensor at this radial slice. As we will not he specifying boundary conditions anywhere, 
we have to parametrize these data in terms of arbitrary functions which will depend on po- 
This parametrization has to be consistent with the hydrodynamic limit, i.e. with the derivative 
expansion. We will see later in this section how these can be achieved by promoting the coeffi- 
cients of the non-normalizable and normalizable modes of the metric fluctuations to arbitrary 
functions of the dimensionless variable poT(x). 

For the boosted black brane solution g^^ = rj^^ and the bd metric can be characterized by 
four parameters: the Hawking temperature Th and the three independent components of the 
unit constant boost four- vector u^, which satisfies u'^rifj^uu" = — 1. The stress tensor is: 



tl^ = IC^^^{v,u + ^u,u,) , where ^ = (2-1) 



and I is the radius of AdS^. 

In the asymptotically slowly varying space-times discussed above, the parameters rj^^, 
and Th can be promoted to arbitrary slowly varying functions of the boundary coordinates x 
131 , p!5| , [50[| . First, we consider a four-dimensional metric h^^{x) such that it reduces to rj^^, far 



in the past or in the future. This metric is also slowly varying such that its curvature scale is 
small compared to the Hawking temperature of the initial or final black brane configuration of 
our space-times. Second, the boost four- vector ■u'^ is promoted to u^{x), an arbitrary slowly 
varying function of the boundary coordinates, which is normalized with respect to the metric 

= -1. (2.2) 

As h^^{x) reduces to ?7^y in the far past or the far future, can be made to coincide with 

the constant black brane boost vector u'^ as well. Also, as we promote m'^ to a function of 
boundary coordinates, without loss of generality, we can use a field definition of such that 
it is the local velocity of energy transport. This can be done by going to the Landau-Lifshitz 
frame locally: 

tl,{x)u''{x) = -e{x)u,{x) , (2.3) 

where e(x) is the local energy density and the index in u^{x) is lowered with the metric h^v{x). 
Now the boundary energy-momentum tensor ( p.l|) is given instead by: 

t(x) = /C^^^^^(V(^) +4«,(a:)«.(x)) , (2.4) 
and the x-dependent temperature T{x) can be defined by: 

e{x) = \lc(7,T{x))\ (2.5) 
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so that T{x) coincides with the Hawking temperature of the black brane in the far past or in 
the far future. 

Notice that the background with the non-constant h^y{x), u^{x) and T{x) will solve the 
5d Einstein equations only to the zeroth order in the derivative expansion. Already at the first 
order in this expansion both the boundary metric and energy momentum tensor will receive 
corrections built from the derivatives of h^^{x), u^{x) and T{x). 

Let us focus on the first derivative correction to the bulk metric. In this paper we will 
consider only the most general tensor perturbations (which will be proportional to the hydro- 
dynamic shear tensor). Here the word "tensor" refers to the transformation properties under 
the spatial rotation group. In our upcoming paper we will study scalar and vector pertur- 
bations showing that the conclusions of this paper will be unaffected. 

As with any other field in the AdS background, there will be two modes for the bulk metric 
perturbation, the normalizable and the non-normalizable one. Let us denote the dimensionless 
coefficients of these modes by Cn and Cnn respectively. 

Let us stress again that throughout this paper these coefficients will be arbitrary, as we 
don't want to specify boundary conditions anywhere. For the sake of completeness, at the end 
of this section we will show how a choice of Cn and Cnn sets specific boundary conditions at a 
given radial slice. 

The non-normalizable mode of the bulk metric does not vanish on the AdS boundary and 
the boundary metric changes for non-zero Cnn [0: 

g^uix) = h,Ax) + + ^(e') ' (2-6) 

where: 

a^,{x) = P^\x)P,''ix){Vxu^ + V^ux){x) - -P^.{x) (V"m„) (x) , (2.7) 

P^iu{x) = h^^{x) + u^{x)u^,{x) , (2.8) 

and the V^^ covariant derivatives are calculated from h^^^x). Also, all the indices are low- 
ered/raised with this metric/it's inverse. Further, here and in all of the upcoming formulae e 
is the derivative expansion parameter. It is the ratio of the typical length scale of variation of 
the boundary fields {h^^{x),T{x),u^^{x)) to the inverse of the late or early time temperature. 

Notice that after considering the full metric perturbation, h^^ will stand for the boundary 
metric only when the non-normalizable mode is switched off. In general it is an useful auxiliary 
variable (along with and T) to obtain the general hydrodynamic metric. 

The boundary energy momentum tensor also receives corrections, but unlike for the metric 
case, the new term is proportional to Ca- 



(7rT(x)) 



4 



C(^) = ^ I '"^V ( V(^) + ^u^ix)u,{x)') - c„ ■ inTix)f a^,ix) + ^(e^) . (2.9) 
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Here (similar to the previous equation) the (nT^x))^ factor in the last term ensures that Cn is 
indeed dimensionless. It is easy to verify that the new tjj^, is both conserved and traceless with 
respect to g^,^, namely: 



g""'" = + Oie') . (2.10) 

In proving this one has to use that a^^ = and u^a^y = 0. Notice that the tracelessness 
condition implies that the boundary fluid is conformal.^ 

We can now present the transport coefficients of the boundary fluid as a function of the 
local temperature and the parameters (cn,Cnn)- To achieve this goal, we have to re- write the 
boundary energy momentum tensor ( |2.9| ) in terms of the proper boundary parameters. First, 
notice that the boundary metric is g^^^x) and not h^u{,x). Second, we have to address the 
normalization of the 4-vector u^{x), which was earlier defined as u^h^'^u^ = —1, while for the 
correctly normalized boundary velocity vector u^^ one needs: 

«^V,.-^^ = -l. (2.11) 

It appears, though, that if we set: 

u^'^ix) =u^{x) + 0{e^), (2.12) 

then the normalization ( p.ll|) will still be satisfied up to the second order in the derivative 
expansion. The physical tensor should also be defined as: 

2 

J 

3 

with: 



^U^) = P''1{x)PK{x) (V^n. + VNa) {x) - \P^,A^) (V^"«„) (x) , (2.13) 



P^ix) = ul{x)u':{x) + gl,{x) . (2.14) 

In the above expressions V'^ is constructed from g^^ix) and all the indices are lowered/raised 
with this metric/it's inverse. Similarly we obtain that: 

cr;, = a,. + 0(e2). (2.15) 

Before proceeding further, let us introduce two new parameters: 

Ci = Cnn and C2 = Cn + ^ . (2-16) 

With the field redefinitions ( p.6|) , ( p. 12 ) and (|2.15 ) the boundary energy-momentum tensor 
(|2.9|) takes the following conformally covariant hydrodynamic form with respect to the metric 

t(^) = /c[^^^^^^(^?J.(x) +4«J(x)^^(x)) - C2 (vrT^x))V^^(x) + 0{e')\ . (2.17) 



: 



■^Recall that the conformal anomaly will show up only at the fourth order in the derivative expansion. 
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Here we set the boundary temperature T^{x) as: 

T^{x)=T{x). (2.18) 
This condition is equivalent to the requirement that the boundary energy density e^{x) will 



have the same dependence on the temperature as in ( p.5|) (this also explains why we have not 
considered any du or dh corrections in the coefficient of the ffist term in (|2.9| )). It follows that: 

^]c(^7rT^{x)y , = ^}c(nT^{x)y , = /C7^(^7rT^(x))^ (2.19) 

where the last two identities follow immediately from ( p. 171) and the thermodynamic relation 
sT = e + p. We also learn from ( p.l7| ) that the shear viscosity is: 

= /Cc2 (7rT^(x))' and !^ = ^ . (2.20) 

We would like to emphasize that in order to arrive at ( p. 171 ) we did not have to use the equations 
of motion. 

To complete this section let us report the full five-dimensional solution that asymptotes to 
the boundary metric ( |2.(j| ) and the boundary energy-momentum tensor ( |2.9| ). 

As we do not yet insist on the regularity condition at the future horizon or at the past 
horizon, and we focus on the renormalization group flow of the energy-momentum tensor, the 
Fefferman-Graham coordinates will be more suitable for our investigations. Further, it is easier 
to work with arbitrary boundary data in the Fefferman-Graham coordinates. Importantly, our 
5c? metric can be translated to the ingoing Eddington-Finkelstein coordinates order by order in 
the derivative expansion, where the change of coordinates is also systematically corrected order 
by order in the expansion [|15|. The details are in Appendix. 

We denote the Fefferman-Graham radial coordinate as p and as before the boundary coordi- 
nates collectively as x. Following |]T3|, we can readily solve the radial flow in these coordinates, 
and we find that up to ffist order in the derivative expansion, the bulk metric takes the following 
form: 

ds^ = — (dp^ + g^^{p, x)dx^'dx''] , (2.21) 

C(p x) 

9f,u{p,x) = A{p,x) h^y{x) +B{p,x) Uf,{x)u^{x) - crf,u{x) + O(e^). 

By definition the functions A, B and C are dimensionless and are given by: 

+ B^f.TWpr^l+'^Ml , (2.22) 



4 / \ \ 1 + (7rT(x)p)V4 
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It is easy to see that ( p.21|) indeed reduces to (|2.6| ) for p = 0, namely at the boundary. On the 



other hand, the horizon is given by (p7rT(x))^ = 4. Notice that the horizon is not a constant 
radial slice, but its radial location varies slowly due to slow variation of the temperature field 
T{x). 

The constraints of Einstein's equations are satisfied by (|2.21|) provided that: 



'"X^ = , (2.23) 



where again the index is raised by g . Besides, had we not imposed ( |2.1U| ) by construction 
of Einstein's equations would have enforced it instead. 



Remarkably, by translating the metric ( p.21| ) to the outgoing/ingoing Eddington-Finkelstein 



coordinates as shown in Appendix, we conclude we have regularity at the past/future horizon 
only when: 

C2 = t\, (2.24) 



b-T— • (2-25) 



so that: 

7]^ _ 1 

We see that it is the physical shear viscosity which alone determines the regularity at the 
past/future horizon. 

Before closing this section let us make a crucial comment regarding the coefficients Ci and 
C2, or alternatively Cn and Cnn- In general, both ci and C2 can be a:-dependent. This dependence, 
however, should appear only through the hydrodynamical variables h^y{x), u^{x) and T{x). In 
our first derivative approximation only T(x) is actually an option. This is so because it takes 
at least one derivative to produce a scalar from h^y{x) and u^{x), and, at the same time, 
Ci and C2 already multiply one derivative terms. This seems to create a contradiction, since 
these coefficients are, by construction, dimensionless and as such cannot depend only on T(x). 
The puzzle is resolved if we assume that the coefficients depend also on a certain value of 
the radial coordinate po- This radial parameter po is a choice of scale, where we specify the 
data (the induced metric and the Brown- York stress tensor) that determines the radial flow. 
The combination poT(x) is dimensionless (and as such preserves conformal invariance at the 
boundary) and thus, still solving the equations of motion, we might have: 

ci = Ci(poT(x)) and C2 = C2{poT{x)) (2.26) 

in ( p.6|) , ( p.9|) and all other formulae in this section. 

As we have already mentioned, regularity flxes only C2 leaving ci completely unconstrained. 
The easiest way to obtain function of pqT{x) is to impose the Dirichlet boundary 

condition on the metric at p = po- In this context the Dirichlet boundary condition implies 
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that the first order fluctuation of the induced metric at the radial slice p = po vanishes. A 
simple inspection of ( |2.21| ) shows that it leads to: 



c.(p„rn.)) = -c.(p„T^.)) log ( '-fy^iif ], (2.27) 



and in conjunction with ( p.24p this determines Ci (T'^(x)po) completely. 



Our approach, however, will be essentially different. We will show that for the fluid/gravity 
correspondence at any hypersurface, both the Dirichlet boundary condition at the cut-off and 
the regularity at the past/future horizon should not be explicitly imposed. In other words, the 
holographic renormalization group flow preserves the fluid/gravity correspondence for arbitrary 
Ci and C2- At any radial hypersurface, with appropriate field redefinitions, we will always get 
unforced hydrodynamic equations exactly as we have got here at the boundary. To arrive at this 
observation, it will suffice to determine only the stress tensor counter-term, without addressing 
the boundary conditions at all. We will see in the next section that with our approximation 
and with some simple additional assumptions this counter-term can be fixed unambiguously. 

With this been said, there is still a possibility that in an implicit way our counter-term 
necessitates a definite boundary condition at a radial hypersurface. We hope to address this 
issue in future research. 



3. Counter-terms and the renormalized energy- momentum tensor 

We recall the traditional prescription for holographic renormalization PSI Bl, 155 



In this prescription, variables like the Brown- York stress tensor are evaluated first at a radial 
cut-off. This is morally dual to the field-theoretic procedure of regularization where we put an 
explicit cut-off to tame ultraviolet divergences. However, these variables, evaluated at the cut- 
off radial slice, diverge at the boundary, and so just like in field theory we need counter-terms. 

In the traditional prescription the counter-terms that renormalize the Brown- York stress 
tensor are derived from a counter-term action living at the cut-off. This is in turn constructed 
by requiring that it cancels the divergences in the on-shell bulk action. However, we should 
also have a well defined variational principle. This is achieved by imposing the Dirichlet bound- 
ary condition at the cut-off. This may require to introduce additional counter-terms like the 
Gibbons-Hawking counter-term which are not related to cancellation of divergences. In order 
to be consistent with the variational principle, the counter-term action required to cancel the 
UV divergences should be a functional of the values of bulk fields at the cut-off, their covariant 
tangential derivatives, but not of their radial derivatives. Further, the traditional prescription 
also requires that the coefficients of all the terms appearing in this counter-term action should 
have the same values for all bulk solutions. 
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As mentioned in Introduction, recent studies indicate that the coefficients of the counter- 
terms may depend on the bulk solution The boundary (counter-term) action is a result 



of integrating bulk degrees of freedom in the region between the boundary and the cut-off, and 
still should be a functional of the values of the bulk fields at the cut-off and their tangential 
derivatives. Nevertheless, it should be obtained in conjunction with the bulk solution in the 
saddle-point action and not just derived from the latter. Thus the boundary action affects the 
boundary condition at the cut-off. 

However, it is not clear to us how to implement such ideas for gravity perturbations in the 
bulk. At the same time, as we do not want to use the Dirichlet boundary condition at the 
cut-off, we also cannot use the traditional prescription of holographic renormalization. Another 
complication stems from the question that, as the metric gets perturbed, whether the correct 
choice of the cut-off hypersurface it still the one which is at a fixed radial location. 

Despite these difficulties, we can make a lot of progress by the following considerations 
without constructing an explicit boundary action: 

• The counter-terms should reduce to the counter-terms of the traditional prescription at the 
boundary. This should be so because the ultraviolet divergences are typically independent 
of the state, for instance the ultraviolet divergences do not depend on the temperature in 
thermal field theory. The coefficients of the counter-terms at the boundary are thus the 
same for all bulk solutions. In other words we have to recover the traditional results at 
the boundary, for instance the Balasubramanian-Krauss energy-momentum tensor. 

• Away from the boundary the coefficients of the counter-terms in our space-times can 
depend on p^T{x) and p'^'V^, with being the covariant derivative constructed from 
the induced metric at the cut-off hypersurface. As in the holographic dual picture, we 
capture not the operator but its expectation value in a state, the coefficients of the 
counter-terms can have a dependence on the state, or equivalently the dual bulk solution. 
Also, p^T{x) ~ T(x)/A^ and p'^V'^ ~ A;/A^ if we make field-theoretic analogies. 

• We introduce all possible counter-terms which depend on the values of the fields at the cut- 
off and their covariant tangential derivatives, even if they are not required for canceling 
ultraviolet divergences. It is plausible that some of the coefficients may vanish at the 
boundary. 

• We require that the renormalized expectation value of the operator takes the same form 
at the cut-off hypersurface after appropriate field redefinitions. Furthermore, these field 
redefinitions should be obtained off-shell, namely without resorting to the equations of 
motion. This is typically true if there is an underlying path integral formalism for renor- 
malization. 
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• Imposing the energy-momentum conservation at the boundary should lead to the analo- 
gous conservation at the cut-off. 

Importantly, our discussion above is general and does not rely on the choice of a renor- 
malization scheme. Instead we try to see what counter-terms one can add on the gravity side 
which conform with expectations on the field theory side. We admit that our discussion does 
not lead to explicit form of the counter-terms in the gauge theory. 

It turns out that the above requirements specify the counter-terms to a large extent. Up 
to first order in the derivative expansion, the only counter-term we can add to the Brown- York 
stress tensor at the cut-off is proportional to: 

/(p^T(x),p<=V^)7^. ~ fip'Tix))^;^, (3.1) 

where we can drop the dependence on p^V^ as the covariant derivative of the metric vanishes. 
Here j^^, is the induced metric at the cut-off p = p^. The next counter-term that we can add 
should involve i?ic[7] or R[y] and hence it contributes two derivative onward. We find this 
counter-term satisfies: 

V'^'^^f {p^Tixj)j;J) = f'{p^T{x))d,T{x) , (3.2) 

where ' denotes differentiation with respect to p^T{x). We find that as long as / {p'^T{x)) is not 
constant the new term violates the energy- momentum conservation. To summarize, / {p'^T[x)) 
should be a constant.^ 

Notice that the physical metric at a given cut-off is not the induced metric itself. In the 
Fefferman-Graham coordinates it is: 

9lv^-j2- %v ■ (3-3) 

It coincides with the boundary metric at p = 0^. Clearly the covariant derivative V 
constructed from the induced metric 7^^ is the same as that constructed from the physical 
metric gl^. 

^Notice that the situation might change if we use other kinds of hypersurfaccs hke p = a/T{x), where a is 
an arbitrary numerical constant. On these hypersurfaccs, we can introduce the counter-term ,f{a)"1^y, where 

now the coefhcient /(a) could be a non-trivial function. As a is a constant we get V'^'' ^/(a)7^j,^ = 0, so we 
do not violate the Ward identity. We require that when a = 0, the coefficient /(a) is fixed by cancellation of 
the ultraviolet divergence in the Brown- York stress tensor. Unfortunately, we find that for such hypersurfaces, 
the renormalized energy-momentum tensor has derivatives of the temperature field, and no field redefinition 
can bring the hydrodynamic energy-momentum tensor to the standard form, unless we use the hydrodynamic 
equations of motion. Thus the renormalized energy-momentum tensor preserves its form on-shell but not off-shell 
as we change a. 

^The multiplying factor above can be replaced by an arbitrary function of pjl vanishing like p^ at 

the boundary. This represents an overall choice of scale and does not affect any dimensionless physical quantity 
like the speed of sound or r}/s at any radial slice. 
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The Brown- York stress tensor is given by: 

yBrown-York ^ _^ (^^^ _ _ ^g ^) 

We emphasize that the Brown- York stress tensor is conserved with respect to the induced metric 
at any hypersurface irrespective of the boundary condition, as this is simply equivalent to the 
pullback of Einstein's equations to the cut-off. Therefore it acts like a bare energy-momentum 
tensor. In the Fefferman- Graham coordinates: 

Up to first order in the derivative expansion, as discussed here, the counter-term for the energy- 
momentum tensor at a constant radial slices is: 

^counter-term ^ _^ ^^3^^^ ^ ^ ^^2) ^ ^ (g g) 

where the coefficient is fixed by requiring the cancellation of the ultraviolet divergence of the 
Brown- York tensor at p = 0. 

The physical renormalized energy-momentum tensor at a given cut-off is: 

V(p, X) = ^ (T^Brown-York ^ ycounter-term^) _ ^3 7) 

The extra multiplying factor of / is determined by the multiplying factor p^/t^ used in 
defining the physical metric in terms of the induced metric in ( ^.3p and is the same as in the 
traditional prescription. 

We find that with the bulk metrics ( p.21|) under consideration, in the Fefferman-Graham 
coordinates the physical energy-momentum tensor up to first order in the derivative expansion 
at the cut-off radial slice is: 



-7(^^) (7rT^x))V(^) + ^(e') 



with: 



«(^) = ./^''L ^ dA^,]- ^d-A{z) , (3.9) 



A{z)-B{z) \A{z 



A{z)-B{z) \A{z) J A{z) 
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where z and and defined by: 

z= (7rr^(x)p)^ and z" = {7iT^{x)p''Y . (3.10) 

However, we still need to find the physical hydrodynamic variables at the cut-off and demon- 
strate that the physical energy-momentum tensor takes the standard hydrodynamic form which 
is covariant with respect to the physical metric at the cut-off. We will do so in the next section. 

4. The RG flow of the fluid 

In the previous section we have found that at an arbitrary radial slice, the physical metric 
g^i,{p,x) is ( |3.3|) and the physical energy- momentum tensor t'^^{p,x) is ( |3.8|) . 

Given that t^^ takes the hydrodynamic form at the boundary, it is certainly not obvious 
that it will take that form at any radial slice. For this purpose we need to determine the 
appropriate hydrodynamic variables at the cut-off, u^'^{x) and T'^{x). 

Following the discussion in Section 2, we define the velocity field with the following 
field definition (analogous to that at the boundary): 

= -e^<7>'^'^ and w'^V^X' = -1- (4.1) 

This definition ensures that u^'^ is the local velocity of energy transport and is normalized with 
respect to the physical metric g^^^ as at the boundary. 

The temperature field at the cut-off T'^ is red-shifted with respect to the temperature at the 
boundary by the factor y/—gfj,iyU'^u'^, which is the square root of the local time-time component^ 
of the physical metric. Therefore: 

^^(P^ x) = + O (e^) . (4.2) 

(a (p^Tb(x)) - B (p^Tb(x)) j 

At the boundary (p^ = 0) we have A = 1 and i? = 0, so clearly T*^ coincides with T^. On the 
other hand, on the horizon {p^T^ = \/2) the red-shift factor diverges and hence so T^. Using the 
explicit form of A and B from ( p. 221) one can express in terms of and p'^. Unfortunately, 
the final expression is not particularly illustrative, as it includes a root of a cubic polynomial, 
and we will not report it here. 

We have to find whether with these field definitions of u^'^ and T'^, the physical energy- 
momentum tensor t^^, at an arbitrary radial slice takes the conformally covariant hydrodynamic 
form with respect to the metric g^^,. 

^Notice that u'^'^d^ is the time-hkc Kihing vector in the un-perturbed black brane background. 



- 15 - 



Using a^j^yU^ = 0, as evident from the definition of a^u, we find that satisfies the first 
equation of ( [4.1|) , but does not satisfy the normahzation condition. This imphes that and 
u'^^ are proportional to each other up to a normahzation factor. We find: 

«^'^(P^ ^) = ^ + O (6^) , (4.3) 



A(p^Tb(x)) - fi(p^Tb(x)) 

where we used (|2.12|) . Non- surprisingly the time-like vector scales exactly as the temperature 
in ( [4.2| ) and the field definition becomes singular at the horizon ^. 

It should he observed that the field redefinitions that take us from the boundary fluid variables 
to the fluid variables at the cut-off are not conformal transformations, as the induced metric at 
the cut-off and the boundary metric are not related by any conformal transformation. We will 
also see that the hydrodynamic energy-momentum tensor at the cut-off is also not traceless 
with respect to the induced metric (up to first order in derivative expansion we can ignore 
considerations of the conformal anomaly, and claim that if the energy-momentum tensor is 
confomal, it should be traceless). 

For the rest of the section we use the convention that whenever the superscript c is shown 
explicitly, the indices are lowered/raised using the metric g^^/iVs inverse. The appropriate 
definition of cr^^ is given by: 

= P\\^)P\\^) (VX + V:<) (x) - ^P;,(x) (V-<) (x) , (4.4) 

with: 

P',A^) = g'^Jx) + u';^ix)u:ix) , (4.5) 

being the projection tensor on the hypersurface orthogonal to with respect to the metric 
g^^, and being the covariant derivative constructed from g^^. Remarkably, using ( [4.3|) , after 
some algebra, we find that: 

A 



;CJ,. + O (6^) . (4.6) 



{A-Bf^ 

The result above follows solely from the definitions of cr^j, and and in proving it we did not 
have to resort to equations of motion. 

It follows from ( |4.3| ) and ( [4. 6] ) that the physical energy-momentum tensor at the cut-off ( |3.(j| ) 
takes the standard covariant hydrodynamic form up to first order in the derivative expansion: 

tl, = }c(a^{p^T%x)) {nT%x)Yg^^,{x) + ^ (p^T^(x)) (7rT^(a;))'n^(x)<(a;) (4.7) 

>3 



-7^(p^T^(x))(7rT^(x))V^,(x) + 0(e2) 



^This indicates that the appropriate four-velocity vector at the horizon is light-like as in the membrane 
paradigm H, 0011,11 . 
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The energy density e*^, the pressure and the shear viscosity rf at the cut-off are: 

= /C (ttT^) ^ - a^) , = /C (^vrT'^) and = K (nT""^ . (4.8) 

Exphcitly: 

.'(/.X) ^ 3K (.r'-(x))' ^'-^i;'!"" = 7^1^^^^ (4.9) 



(1 -^74) (1 + ^74) 

(1 + ^74)3/2 



where z and are defined as in ( p. 10 ). To derive these formulae we first plugged ( |2.21| ), 
and ( |4.6| ) into ( [4.7| ), then compared it with ( p.8| ) and ( |3.9| ) in order to express e'^, and r^'^ in 
terms of A{z), B{z) and C{z). Finally, we substituted the explicit solution of the equations of 
motion from ( |2.22| ) to arrive at the final result. It is important to emphasize here that T^{x) 
on the right hand side of all three equations in ( [4.9|) should be seen as a function of T'^{x) 
and p*^ (both in the overall factors and in the definition of z^), since we want e^, ]f and Vj'^ to 
depend only on the cut-off parameters p'^ and T'^(x). The reason for this implicit dependence 
on T'^{x) is the complicated inverse of the function given in (|4.2| ) (see the discussion below this 
equation). 

It is worth to notice that both and t]^ are finite at the horizon {z'^ = 4), while the local 
pressure diverges (see |1^, [1^ for the earlier discussions). 

From ( [4. 9] ) and the thermodynamic relation e'^ + = T'^s'^ we arrive at our final result: 



Comparing with ( |2.2(J| ), we learn that: 



!L.£l. (4.10) 

TT 



^ = ^ (4.11) 



Thus, as advertised in Introduction, the physical ratio r]/s obtained after appropriate field 
redefinitions, does not flow radially irrespective of the explicit choice of data at po that com- 
pletely determines the radial flow (or equivalently the explicit choice of boundary conditions at 
the cut-off and at the horizon). Further, requiring regularity in the past/future horizon ( p. 241) 
implies that: 

1 

^ = T^- (4-12) 

Air 

Thus the above is also true irrespective of the explicit boundary condition at the cut-off. This 
result has been foreshadowed by the study of scale dependence of the holographic response 
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function \El 



|. It has also been derived in p6| in the context of holographic non-relativistic 
hydrodynamics. Interestingly, again as in in the context of the holographic non-relativistic 
hydrodynamics, the bulk viscosity also vanishes at all scales. 
Additionally, we find that: 

~ (l + zc/4)3/2 ' (4.13) 

therefore the entropy density decreases monotonically as move away from the boundary to the 
horizon. This is expected from the general renormalization group flow arguments. The total 
entropy and not the entropy density is specified by the area of the horizon. In the case of planar 
black holes the area of horizon is infinite, so the entropy density should scale as inverse of the 
3d area. This is exactly what we find. 

To conclude, let us present the result for speed of sound at the cut-off, which follows directly 
from (^191): 

t^sound 3 l_2c/4 " ^ ' ^ 

We see that Csound blows up at the horizon and is equal to 1/a/3 on the boundary in accordance 
with [0, |1^, |23|, |2^, ^ . The divergence of the speed of sound at the horizon (2;^ = 4) is most 
likely related to the incompressibility (V ■ m = 0) of the fluid on the horizon 



5. Discussion 



We have found that the renormalization group flow of energy-momentum tensor at a radial slice 
of space-time as given by Einstein's equations preserves the hydrodynamic form of the energy- 
momentum tensor up to first order in the derivative expansion for arbitrary slowly varying 
boundary data which equilibrate in the far past or far future. This result is kinematic, in the 
sense that to show this we only need to solve for the radial fiow of the boundary data, without 
imposing the constraints implied by Einstein's equation on it. Also, this result is independent 
of the explicit choice of boundary conditions for the radial flow. Our results thus shed some 
light on holographic renormalization in general. 

Much of our procedure can be generalized to higher derivative corrections to the Einstein- 
Hilbert action which are local functionals of the Riemann curvature and its covariant derivatives. 
However, we feel that we need to understand the renormalization group flow given by Einstein's 
equations better flrst. We mention a few open questions in this context: 

1. It is important to flgure out whether the hydrodynamic form of the renormalized energy- 
momentum tensor is preserved off-shell by the radial flow of the boundary at second and 
higher orders in the derivative expansion as well. We expect this will help to address 
the general question of the choice of counter-terms and hypersurfaces involved in deflning 
holographic renormalization. 



- 18 - 



2. What are the degrees of freedom at the horizon? In particular, we have to understand 
if the hydrodynamics at the horizon is purely incompressible Navier-Stokes. To address 
this problem, one needs to study the higher derivative corrections to the bulk metric. In 
principle, the answer can also depend on the renormalization scheme. 



3. What is the renormalization group flow of the non-hydrodynamic part of the boundary 
energy-momentum tensor? Non-hydrodynamic data at the boundary in the context of 



AdS/CFT has been investigated in [^8|, 0. The study of the renormalization group flow 
of such non-hydrodynamic data will help to answer the question if non-hydrodynamic 
degrees of freedom are relevant at a regular future horizon. 



The study of RG flow of the energy-momentum tensor in the hydrodynamic regime itself 
may throw light on the holographic emergence of gravity as argued in |^3|, ^ ^ ^ . In 
particular, we know if that if the transport coefficients of the holographic fluid are not 
right, we end up with naked singularities. If naked singularities are unobservable, we 
need to understand this in the context of holographic emergence of gravity. It may be 
possible that absence of naked singularities comes from appropriate physical restriction 
on the renormalization group flow of the boundary data. 
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A. Translation to Eddington-Finkelstein coordinates 

The metric can be translated to the ingoing/ outgoing Eddington-Finkelstein coordinates order 
by order in the derivative expansion. The coordinate transformations also receive corrections in 
the derivative expansion which can be systematically computed. We will denote the Eddington- 
Finkelstein radial coordinate as r and the boundary coordinates of this coordinate system 



collectively as y. Following |15|, we flnd that up to flrst order in the derivative expansion, in 
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the ingoing Eddington-Finkelstein coordinates the metric (|2.21|) is: 



Gf,y{y,r) = r'^hf.^iy) + (^^y) u^{y)u^{y) 



nr + 7iT{y)f{r^ + {nT{y)y) 
27cT{y) \ H 



— 2 arctan 



(A.l) 



.T^T{y) 



~r I {u{y) ■ V){Uf,{y)u^{y)) - -u^{y)u^{y){V ■ u){y) 



+Cl 



Tiy) 



+0 (e^) 



r'^cr^^{y) + C2 



T{y)\ 1\ 



A 7rT{y) 



log 1 



Wy)y 



where: 



Cl 



Ci(poT(x)) + 0(e) 



C2 



C2(poT(x)) + 0(e) 



(A.2) 



parametrizing the choice of data at the Eddington-Finkelstein radial slice ro just as at the 
Fefferman- Graham radial slice po in Section 2. 

In the above coordinates the horizon is at r = 7rT(?/), where G^uU^{y)u'^{y) vanishes. Also 
the boundary data {u^{y),T{y), h^^{y)) are the same functions of new boundary coordinates 
y. The boundary is located at r = cxd where the new boundary coordinates y coincide with the 
Fefferman- Graham boundary coordinates x, though they are nontrivially related in the bulk. 
Except for the penultimate line in ( |A.1| ), this metric agrees with that in where arbitrary 
perturbations of the induced metric and the Brown- York tensor have not been considered as 
here. 

We also see from the explicit form of the metric in the ingoing Eddington-Finkelstein 
coordinates that the future horizon is regular when C2 = 1/4, i.e. when rj^/s'' = 1/{Att). The 
metric in outgoing Eddington-Finkelstein coordinates can be obtained by reversing the sign 
of the time coordinate and also the time component of u^{y). We readily observe that for 
regularity of the past horizon, we need yf / = —1/{Att). The reversal of the sign simply 
follows from the time reversal invariance of Einstein's equations. 
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